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We have performed /iSR studies on single crystals of SrCu2(B03)2, a quasi-two-dimensional spin 
system with a spin singlet ground state. We observe two different muon sites which we associate 
with muons located adjacent to the two inequivalent O sites. One site, presumed to be located in 
the Cu-O-Cu superexchange path, exhibits a large increase in the frequency shift with decreasing 
temperature which is unaffected by the singlet formation, indicating that the muon has locally 
broken the singlet bond. We have also performed /iSR on single crystals of SrMgo. 05 Cui. 95(603)2, 
Sro.96Lao.o4Cu2(B03)2, and Sro.9sNao.o5Cu2(B03)2. We have found that the frequency shifts of 
these doped samples are equivalent and contain three branches at low temperatures. Two of these 
branches map on to the branches observed in the pure sample reasonably well and so we attribute 
the third branch to effect of the dopants. Specifically, this third branch represents the case when 
the muon sits at a site in the superexchange path lacking a corresponding singlet, due to it already 
being broken as a result of doping. This then leads to the conclusion that singlets are broken in this 
system when it is doped both in and out of the CUBO3 planes. 

PACS numbers: 76.75.-K, 75.50.-y 



I. INTRODUCTION 



Low-dimensional spin systems with spin-singlet ground 
states have attracted considerable study over the past 
15 years. These systems are interesting for a number 
of reasons, including possible relevance to high tem- 
perature superconductivity and the existence of novel 
ground states. Most of the systems studied to date 
have been limited to quasi-one-dimension. These include 
spin ladder systems (e.g. SrCu203[l|), spin-Peierls sys- 
tems, (e.g. CuGeO3 0]) and Haldane gap systems (e.g. 
Y2BaNiO5 0]). These materials all have a spin gap in the 
energy spectrum, which is associated with the spin sin- 
glet ground states. Theoretical studies of spin ladders[4( 
have predicted d-wave superconductivity at low concen- 
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trations of hole doping in these systems. In fact, the 
mixed ladder-chain system Sri4_ a; Ca a ;Cu2404i has been 
found to be a superconductor under pressure, specifically 
in the cases of x=13.6[5j and x=11.5,[6[ with the super- 
conductivity thought to originate in the ladder layers. 
High- T c cuprates have d-wave superconducting states 
and a partial (or pseudo) spin gap, which provides fur- 
ther evidence for strong similiarities between doped lad- 
ders and doped Cu02 planes and motivates study of the 
effects of perturbing spin gap systems. 

Two dimensional systems with spin singlet ground 
states are comparatively rare and often involve lat- 
tices where the active spins have a small number of 
nearest neighbours, such as CaV40g,[3| where the de- 
pleted square lattice consists of plaquettes of four spins, 
coupled to neighbouring plaquettes through one ver- 
tex. SrCu2 (603)2 is a second example of a quasi-two- 
dimensional spin system with a spin-singlet ground state. 
Q The unit cell is tetragonal and the structure is char- 
acterized by layers consisting of Cu 2+ , 2 ~ and B 3+ 
running perpendicular to the c-axis. These layers are 
separated from each other by planes composed of Sr 2+ 
ions. All Cu 2+ sites have localized spin S = 1/2 mo- 
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ments and one nearest neighbour with which they form 
dimeric units, while each of these units is surrounded 
by four orthogonal dimers. This configuration of spins 
corresponds to the two-dimensional model that Shastry 
and Sutherland first introduced in 1981. |9| This model 
is topologically equivalent to a two-dimensional Heisen- 
berg model with nearest-neighbour coupling J and next- 
nearest- neighbour coupling J'. With decreasing J'/ J, a 
phase transition from a four-sublattice Neel state to a 
gapped spin singlet state was shown to take place at 
a critical value of J'/J ~ 0.7. The ratio of J'/J in 
SrCu2 (B03)2 is about 0.68, a little bit smaller than 
the critical value. ^fj| High resolution inelastic neutron 
measurements [ll| have identified one, two, and three 
triplet excitations. These measurements clearly show the 
appearance of the energy gap in the spectrum of excita- 
tions with decreasing temperature, and from them the 
energy gap is evaluated to be ~ 34 K. The one and two 
triplet excitations both display a very large drop-off in in- 
elastic intensity with increasing temperature. This tem- 
perature dependence is generally not seen in a system 
undergoing a phase transition, and no evidence of sym- 
metry breaking has been associated with the appearance 
of the spin singlet ground state. Another interesting fea- 
ture of this system is the magnetization plateaus that 
are observed at 1/8, 1/4,Q and l/3[l2j of the saturation 
magnetization in high magnetic fields (> 20 T). This is 
the first 2D system where these plateaus have been ob- 
served. 

While superconductivity has been predicted and ob- 
served in some quasi-one-dimensional systems with a 
spin singlet ground state, it has also been predicted in 
the doped Shastry-Sutherland model. [L3L ITS Il5| Liu et. 
al have recently investigated a series of polycrystallinc 
SrCu2 (1303)2 samples, with dopants added both infill] 
and out of the CUBO3 planes. They found that the 
spin gap behaviour was strongly affected by doping, but 
did not find any evidence for superconductivity. The re- 
sistivity of the pure sample at room temperature was 
found to be 41.73 Mf2m,[l(| and the resistivity of all 
doped samples was found to be only one or two orders 
of magnitude smaller. This implies that all their samples 
were insulators and that the carriers introduced as a re- 
sult of doping were not very mobile, at least at room tem- 
perature. However, the theoretical predictions discussed 
above make it important to better understand Shastry- 
Suthcrland systems such as SrCu2(B03)2 and their re- 
sponse to perturbations. 

Muon spin relaxation measurements of a spin sin- 
glet state give a signal characteristic of a non-magnetic 
state. Zero field measurements [l?) of the spin-Peierls sys- 
tem CuGe03 gave essentially a temperature-independent 
signal, with weak relaxation being ascribed to isolated 
defects associated with impurities. In the case of 
SrCu 2 (B0 3 )2, Fukaya et al.^ reported zero field (ZF) 
and longitudinal field (LF) fj.SH measurements of ceramic 
specimens. In contrast to their expectations of seeing no 
spin relaxation associated with this state, the authors re- 



ported the onset of a large, roughly constant relaxation 
for temperatures of 3 K and below. They attributed this 
relaxation to dilute defects or to /iSR detecting weak in- 
trinsic magnetism that other probes could not see. The 
authors also reported that dynamic fluctuations persisted 
in the system down to base temperature (20 mK) and 
that these fluctuations likely have a quantum origin. 

Lappas et al. also investigated ceramic samples of 
SrCu2 (1303)2 using the /iSR technique. fi~9j The authors 
reported the same anomalous relaxation in ZF as ob- 
served by Fukaya et al. using a different batch of sam- 
ples. They attributed this relaxation to an unusual spin 
freezing process. The authors went on to perform 0.6 
T frequency shift measurements on their samples. They 
concluded that implanted muons may liberate spin den- 
sity in the spin gap regime that spin-freezes at very low 
temperatures. 

II. EXPERIMENTAL DETAILS 

We grew single crystals of SrCu2 (603)2, 
SrMgo.o5Cui.9 5 (B0 3 )2, Sr .96Lao.o4Cu2(B0 3 )2, and 
Srn.3RNan.nflCu2(BQ3)2 in floating zone i mag e furnaces 
at Tokyo [2(| and McMaster University [21| following 
slightly different procedures. Neutron scattering mea- 
surements have been reported by Haravifard et al.[22j 
on the system SrMgo.05Cu1.95 (603)2- Aside from the 
original crystal growth article [2l|, this seems to be 
the only paper describing doped single crystals of this 
system as they have proven to be quite difficult to grow. 

DC susceptibility measurements of a McMaster crystal 
of SrCu2 (603)2 in a field of 5 T parallel to the c-axis 
are shown in Fig. [TJ results for the Tokyo crystals are 
essentially identical. There is a sharp decrease in the 
susceptibility below 15 K which provides evidence for the 
spin singlet ground state. A small upturn is seen in the 
susceptibility curve below 4 K, which likely arises due 
to magnetic impurities and/or defects of Cu 2+ ions in 
SrCu2 (603)2- The susceptibility data below 7.5 K was 
fit to the following functional form: 

X = C/T + Be-^ T + X0 (1) 

The first Curie term is due to the contribution of mag- 
netic impurities, the second term is to account for the 
spin gap, and the third temperature-independent term 
accounts for Van-Vleck paramagnetism and core diamag- 
netism. The spin gap A was found to be ~ 20.2(7) K in 
our field of 5 T, which is in line with other reported es- 
timates. If we ascribe the Curie term to free Cu spins, 
then this concentration would be only 0.12(4) %. 

We performed 5 T transverse field (TF) ^SR measure- 
ments on the pure and doped SrCu2(60s)2 samples to in- 
vestigate the presumed non-magnetic spin singlet ground 
states in these systems. These measurements were con- 
ducted on the M15 surface muon channel at TRIUMF in 
Vancouver, Canada, using a helium gas flow cryostat in 
the temperature range 2 K < T < 150 K. The samples 
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FIG. 1: DC volume susceptibility of single crystal 
SrCu2 (B03)2, measured in an applied magnetic field 
H = 5 T II c. 



J 







678 

Frequency (MHz 



FIG. 2: (Color online) Fast Fourier transform of muon preces- 
sion signal for SrCu2(BOs)2 measured in an applied magnetic 
field H = 5 T || c at 20 K. The two peaks correspond to two 
magnetically inequivalent muon sites. 



were mounted in a low-background spectrometer with the 
c-axis parallel to the incoming muon momentum (and the 
external field). In this high field experiment, the mea- 
surements of the muon precession frequency signal in our 
samples and a reference material (CaCOs) were taken si- 
multaneously to allow for a precise determination of the 
frequency shift. 

In addition, we performed ZF-/iSR and low TF-/iSR 
measurements (H < 1 T) on the pure sample. After 
observing the anomalous relaxation previously seen by 
Fukaya et q/.[l8| in the powder sample, we performed 
LF-/iSR on our single crystal to characterize the nature 
of this relaxation. For these low temperature measure- 
ments (between 20 mK and 10 K), we used an Oxford 
Instruments dilution refrigerator on the M15 beam line. 
The sample was mounted on a silver sample holder, with 
the c-axis parallel to the incoming muon beam, as in the 
case of the high TF measurements. 



III. DISCUSSION AND ANALYSIS 

Implanted positive muons reside at sites which are min- 
ima in the electrostatic potential. Generally in oxides, 
muons are hydrogen-bonded to O ions, roug hly 1 A 
away. Fig. [3 shows a fast Fourier transform of the trans- 
verse field muon precession signal measured in an applied 
magnetic field H = 5T||catT = 20K. Two peaks 
are clearly visible which correspond to two magnetically- 
inequivalent muon sites in SrCu2(BOa)2; we associate 
these with the two magnetically-inequivalent oxygen sites 
in the system. To take these two muon sites into account, 
this 5 T TF data was fit to the sum of two exponentially- 
damped cosine functions. Our exact fitting function is 
given by: 



P(t) = Ax exp~ Alt cos(LUit+4>i)+A 2 exp 



" Aat cos(w 2 i+<fe) 
(2) 

where oj\ and 102 are the average muon precession fre- 
quencies, Ai and A2 are the muon relaxation rates for 
the two sites, and A\ and A2 represent the relative oc- 
cupancy of the two unique muon sites. Ai and Ai are 
both temperature independent and were found to have 
a ratio of roughly 3:2. The ratio of the two inequivalent 
O sites in the unit cell is 2:1, with fewer O sites lying in 
the Cu-O-Cu superexchange path. However, the prob- 
ability of occupying a specific site depends not just on 
the relative number of such sites, but also on the exact 
energies of those sites. Generally, there should only be 
one true stable muon site, but a metastable site could be 
long-lived on the time range of our experiment, leading 
to the site-occupation ratio that we observe here. 

The frequency shift is the fractional difference between 
the applied field and the local field at the muon site, and 
it is proportional to the local spin susceptibility. We see 
that for one of the muon sites the frequency shift roughly 
follows a Curie- Weiss temperature dependence for T > 4 
K (represented by the crosses in Fig.[3ji), with the precise 
value of the Curie constant depending on the tempera- 
ture range being fit. This is in contrast to 11 B-NMR 
frequency shift measurements. [23| In that case, the au- 
thors applied a field of 7 T perpendicular to the [1-10] 
direction and observed a decrease in the frequency shift 
below the singlet formation temperature of 15 K. The 
lack of a similar decrease in our /i + frequency shift for 
this site indicates the presence of a larger local field than 
expected in the singlet regime. The fundamental differ- 
ence between the two experiments is that in our case 
there is a muon present in the system, which means the 
larger local field is likely the result of a muon perturba- 
tion effect. More specifically, it seems that the muon has 
locally broken a spin singlet bond and created at least 
two quasi- free spins. These spins are then polarized by 
the applied field which leads to our observed frequency 
shift. As this appears to be a strong effect, it then seems 
reasonable to associate this muon site with the oxygen 
site in the Cu-O-Cu superexchange path. 

The shift at the other muon site is considerably smaller 
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FIG. 3: Frequency shifts vs. temperature measured in an applied field of 5 T for: (a) SrCu2(B03)2, (b) Sro.96Lao.o4Cu2(BC>3)2, 
(c) SrMgo.osCui. 95(1303)2, and (d) Sro.95Nao.osCu2(B03)2. The insets show the temp, dependence of the frequency shift for 
the muon site outside of the superexchange path. 



and negative (see inset of Fig. [3^) . In this case, the fre- 
quency shift follows the temperature-dependence of the 
bulk susceptibility reasonably well, indicating the lack of 
any strong muon perturbation. For this reason, this site 
likely corresponds to the O 2- ion out of the Cu-O-Cu 
superexchange path. 

Muon perturbation effects have been observed in other 
materials. In the system PrNi5, it has been shown that an 
implanted muon locally alters the crystal field levels. 24 
Muon-induced singlet-breaking, similar to what is being 
described in this work, has been observed in the spin-1/2 
ladder system KCuCi3.[25| An implanted muon has also 
been reported to perturb a link in the AF chain of the 
quasi 1-D systems KCuF3[26[ and dichlorobis (pyridine) 
copper (II) [27| , where the observed local spin susceptibil- 
ity agreed well with theoretical calculations. [28l |29| Our 
measurement of the local susceptibility at two different 
sites in SrCu2 (1303)2 should serve to stimulate similiar 
calculations for the two-dimensional Shastry-Sutherland 
system. 

Motivated by the discovery of the two muon sites in 
SrCu2(B03)2, we fit the ZF-yuSR data to the sum of two 
terms, each taken to be the product of an exponential and 
a static Gaussian Kubo-Toyabe relaxation function. The 
Gaussian Kubo-Toyabe part represents relaxation due to 
nuclear moments that are oriented randomly in ZF. This 
is generally a temperature-independent contribution to 
the relaxation and this was assumed in this case by fixing 
the Gaussian Kubo-Toyabe relaxation rates. The ratio of 
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FIG. 4: Exponential relaxation rates vs. temperature for zero 
applied field in SrCu2(BOs)2. 



the asymmetries of the two sites was fixed to the value 
obtained when fitting the high TF data. Fig. 4 shows the 
exponential relaxation of one site is roughly independent 
of temperature for T > 4 K, then rises sharply at 4 K and 
saturates for temperatures of 3 K and below. The other 
site has a very small, almost temperature-independent 
relaxation. The large increase and saturation of the ex- 
ponential relaxation for the one muon site at low temper- 
atures may imply slowing down or freezing of the spins 
liberated by the muon perturbation effect as previously 
suggested by Lappas et aL[lj| 



5 



TABLE I: Estimated spin gaps, Curie constants, and fraction of free spin 1/2-Cu's by fitting susceptibility data with H = 1 
kG || ab plane. 



Sample 


Spin Gap (K) 


Curie Constant (K) 


Estimated fraction of free Cu spins (%) 


SrCu 2 (B0 3 ) 2 
SrMgo.osCui.95(B0 3 )2 
Sro.96Lao.o4Cu2(BC>3)2 
Sr .95Nao.o5Cu2(B0 3 )2 


27.8(8) 
20.4(6) 
23.2(6) 
24.4(5) 


3.4xl0 _a (8) 
5.5xl0 _4 (l) 
3.09xl0~ 4 (l) 
1.67xl0" 4 (l) 


0.0585(1) 
0.95(2) 
0.53(2) 
0.29(2) 
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FIG. 5: (Color online) LF-/1SR line shapes for selected applied 
fields at 20 mK. There is significant relaxation even in the case 
of 2 kG. The lines are a guide for the eye. 
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FIG. 6: (Color online) Fast Fourier transform of precession 
signal for SrMgo.o5Cui.9s(B03) 2 measured in an applied mag- 
netic field H = 5 T || c at 20 K. There are now three peaks 
instead of two; the third peak is due to the Mg 2+ ions. 



To see if this was a suitable description, it was neces- 
sary to examine longitudinal field (LF) fj.SH data taken 
below the saturation temperature, specifically at 20 mK 
in this case (Fig. [5]). Values of ~ 5 G and 80 G were 
estimated for the local fields at the two muon sites from 
the ZF time spectrum at T = 20 mK using the rela- 
tion B ~ (7 P = 135.5 MHz/T) with A represent- 
ing the muon relaxation rate as before. Assuming these 
local fields are static, a comparable LF (up to one or- 
der of magnitude larger) should be enough to completely 
decouple them. However, we are still seeing significant 
relaxation even in the case of LF = 2 kG. This implies 
that the spins are dynamic in the low temperature regime 
and their behaviour cannot be described by spin freezing. 
These spin fluctuations are fast for T > 4 K and then slow 



down substantially below T~4K. The saturation of the 
relaxation at the low temperatures suggests that these 
fluctuations are not thermal but instead have a quantum 
origin. This is in agreement with what Fukaya et al. re- 
ported in the ceramic samples. [18| If the LF spectra are 
examined closely at early times, it becomes clear that 
they have a distinct Gaussian nature. In general, dy- 
namic spin fluctuations cannot be described by a Gaus- 
sian lineshape. This appears to be one of the rare cases 
of undecouplable Gaussian behaviour in /iSR, which is 
still poorly understood. Such behaviour has been pre- 
viously observed in frustrated and/or low-dimensional 
systems. [13, HI| 

Fig. [5] shows a fast Fourier transform of the trans- 
verse field muon precession signal measured in an ap- 
plied magnetic field H = 5 T [| c at T = 20 K for 
SrMgo.o5Cui. 95 (B03)2. Similar FFTs were obtained for 
the La and Na-doped samples. There are clearly three 
frequency peaks in this case as opposed to the two that 
were seen in the FFT of the pure sample. For this rea- 
son, the data was fit to the same functional form as given 
by Eq. (2), but three exponentially-damped cosine terms 
were included instead of two. The frequency shifts calcu- 
lated from this fit are given in Fig. [3}>d; the shifts for all 
three doped samples examined appear to be essentially 
equivalent. The upper and lower branches are remark- 
ably similar to the two branches in the pure case, but the 
middle branch is new and so it seems is associated with 
the dopants. 

In one situation, we are substituting non-magnetic Mg 
in for the Cu sites. It is well-known that this leads to 
the direct breaking of spin singlets, as some Cu spins no 
longer have a nearest neighbour Cu with which to form 
a dimer. In the other two cases, we are substituting La 
or Na in for the Sr sites. These are both forms of charge 
doping, with each La introducing one electron into our 
system and each Na introducing one hole. Each elec- 
tron or hole introduced by doping is expected to break 
a singlet pair and make a localized spin- 1/2. Resistivity 
measurements performed by Liu et al. [101 ] on La and Na- 
doped polycrystalline samples verify the localized nature 
of these carriers. The breaking of singlets by the different 
dopants is supported by DC susceptibility measurements 
performed on our single crystals with an applied field H = 
1 kG |j ab plane. [2l[ All three doped samples had much 
larger Curie tails than the pure sample, indicating the 
presence of many more free spin- 1/2 Cu's in the singlet 
regime. We have calculated the spin gap, Curie constant 
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and the resulting fraction of free spins by fitting this data 
to Eq. (1); the results are shown in Table I. In all the 
doped cases, the spin gap has been supressed. Also, the 
fraction of free spins is at least five times as great as in 
the pure case, but still lower than what we expect con- 
sidering the concentration of our starting materials. This 
is likely due to difficulties involved in doping this system 
and suggests that our doping concentrations should be 
viewed as nominal values. 

Our /xSR results complement the susceptibility mea- 
surements by providing microscopic evidence that sin- 
glets are broken in this system when it is doped both 
in and out of the C11BO3 planes. It appears there are 
three distinct possibilities for muons in the doped sam- 
ples. The muon can either stop at the sites in or out 
of the superexchange path as before. If the muon stops 
at the site in the superexchange path, it will break the 
singlet and create two quasi-free spins in close proximity. 
However, in some cases it seems the singlet will already 
be broken due to the dopants and then the muon will 
have no singlet to break. In this case, the muon will 
only be in close proximity to one quasi-free spin and so a 
smaller frequency shift will be measured with decreasing 
temperature, giving us the middle branch. In fact, Fig- 
ure [3)>d shows that this middle branch gives about half 
the shift of the upper branch. 



We also observe a large relaxation increase for the 
site in the superexchange path in the singlet regime, 
which is in contrast to the weak, practically temperature- 
independent relaxation expected in such a state. Anal- 
ysis of LF-^SR data has shown that this relaxation is 
likely due to slow but dynamic fluctuations of the liber- 
ated spins that persist down to our base temperature of 
20 mK. 

Finally, all three doped samples studied have remark- 
ably similar frequency shifts. In each case, there are now 
three branches at low temperatures instead of two as ob- 
served in the pure sample. The upper and lower branches 
map on to the two branches in the pure sample very well, 
and so the third branch is attributed to the dopants. 
Specifically, this third branch represents the case when 
the muon sits at a site in the superexchange path lacking 
a corresponding singlet due to it already being broken 
as a result of doping. Our results demonstrate that elec- 
tron or hole doping of SrCu 2 (1303)2 by substituting on 
the Sr site gives the same result as substitution of Cu 
by non-magnetic Mg: namely, the breaking of a localized 
singlet and the liberation of a single Cu spin. We find 
that the number of free spins is significantly less than ex- 
pected from the nominal dopant concentration and that 
reaching doping levels sufficient for metallic (and super- 
conducting) behavior may be difficult in this system. 



IV. CONCLUSION 

We determined there are two different muon sites in 
SrCu2(B03)2, which we associate with muons located ad- 
jacent to the two inequivalent O sites in the system. One 
site exhibits a large increase in the frequency shift with 
decreasing temperature that is unaffected by singlet for- 
mation, indicating that the muon sits in the Cu-O-Cu 
superexchange path and has broken the corresponding 
dimer. The frequency shift of the other site is consider- 
ably smaller and negative. It roughly scales with the DC 
bulk susceptibility, indicating the lack of a perturbation 
effect in this case. 
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